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• Premise of the study: Bathysa australis is a common subcanopy tree from the Atlantic Forest that is pollinated by bees and 
wasps and produces autochoric seeds. This species exhibits great phenotypic plasticity along the elevational gradient of Serra 
do Mar in southeastern Brazil. We expect to assess the genetic diversity and gene flow between populations of this species 
along the elevational gradient. 

• Methods and Results: We developed a micro satellite-enriched genomic library for B. australis, and 10 microsatellite loci were 
successfully amplified, ranging from one to 13 alleles per locus. The observed and expected heterozygosities ranged from 0.333 
to 0.900 (average: 0.629) and 0.564 to 0.900 (average: 0.742), respectively. 

• Conclusions: These are the first microsatellite markers developed for the genus Bathysa and may be useful in other species of 
the Condamineeae tribe. These primers will be an important tool for studies of population ecology and conservation genetics. 

Key words: Atlantic Forest; Bathysa australis; conservation genetics; medicinal plant; polymorphism; population ecology. 



Bathysa australis (A. St.-Hil.) Hook. f. ex K. Schum. (Rubi- 
aceae) is a subcanopy tree that is widespread along the elevational 
gradient (100-1000 m a.s.l.) of the Atlantic Forest of Serra do Mar 
in Sao Paulo State, Brazil. This species is a common plant in parts 
of the Atlantic Forest (e.g., Ramos et al., 201 1) and has an impor- 
tant role in ecosystem functioning, e.g., providing a nectar source 
to a variety of insects (Andrich, 2008). Furthermore, its bark is 
used in folk medicine (Germano-Filho, 1999), which indicates its 
social value in addition to its ecological value. Bathysa australis 
displays great phenotypic plasticity in leaf size and color along its 
elevational gradient, is pollinated mainly by bees and wasps 
(Andrich, 2008), and presents autochoric seed dispersal (Pedroni, 
2001). Because we believe the elevational gradient might function 
as a barrier for B. australis gene flow, the investigation of the spa- 
tial distribution of the dispersal, pollination, and genetic diversity 
of this plant could generate important information regarding its 
population biology. Above all, the Brazilian Atlantic Forest is sig- 
nificantly threatened (Myers et al., 2000), and the microsatellites 
developed in this study should serve as tools to evaluate the im- 
pacts and define conservation strategies. 
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METHODS AND RESULTS 

We extracted genomic DNA from leaf tissue samples using the DNeasy Plant 
Mini Kit (QIAGEN, Valencia, California, USA). A microsatellite-enriched library 
was then developed following Billotte et al. (1999). The DNA samples were di- 
gested using the Rsal restriction enzyme (Invitrogen, Carlsbad, California, USA) 
for 3 h at 37°C, and the resulting fragments were ligated to Rsal adapters for 2 h at 
20°C. The fragments containing microsatellites were selected by hybridization with 
(CT)8- and (GT)8-biotinylated oligonucleotides, followed by capture with Strepta- 
vidin MagneSphere Paramagnetic Particles (Promega Corporation, Fitchburg, 
Wisconsin, USA). The selected DNA fragments were PCR-amplified in a final 
volume of 100 jiL containing 20 jiL of selected fragments, Ix PCR buffer, 1.5 mM 
MgCl2, 200 |LiM dNTPs, 4 pmol of primer Rsa21, and 2.5 U of Tag DNA poly- 
merase. A PTC- 100 thermal cycler (MJ Research, Waltham, Massachusetts, USA) 
was used with the following program: 95°C for 1 min, followed by 25 cycles of 
denaturation at 94°C for 40 s, 60°C for 1 min, extension of 72°C for 2 min, and a 
final extension of 72°C for 5 min. The amplification products were cloned into the 
pGEM-T Easy Vector (Promega Corporation). Plasmids were transformed into 
Escherichia coli XLl-Blue competent cells, and positive clones were selected 
using the p-galactosidase gene and grown overnight in an HM/FM medium with 
ampicillin. A total of 96 positive clones were bidirectionally sequenced using an 
automated ABI PRISM 377 sequencer (Applied Biosystems, Foster City, Califor- 
nia, USA) with T7 and SP6 primers and the BigDye Terminator version 3.1 Cycle 
Sequencing Kit (Applied Biosystems). The sequences were assembled and edited 
using SeqMan Pro (DNAStar Inc., Madison, Wisconsin, USA). The repetitive 
regions were identified using the Simple Sequence Repeat Identification Tool 
(Temnykh et al., 2001), and 30 primer pairs were designed using WebSat 
(Martins et al., 2009). Ten primer pairs ampHfied microsatellite regions and 
were selected for screening (Table 1). The remaining loci were discarded due to 
amplification failures or nonspecific amplification patterns. The forward primer 
for each pair was labeled with fluorochromes (HEX and TET). 

PCR ampHfications were performed in a 15-jiL volume containing 15 ng 
DNA, Ix PCR buffer, 0.15 mM each dNTP, 0.8 mM each primer, 0.04% bovine 
serum albumin (BSA), 1.5 mM MgCl2, and 1 U Taq DNA polymerase. A PTC- 
100 thermal cycler (MJ Research) was used with the following program: 96°C for 
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Table 1. Characteristics of 10 microsatellite markers developed for Bathysa australis. 

Locus Repeat motif Primer sequences (5'-30 (°C) Allele size range (bp) GenBank accession no. 



BA02 


(CT)8 


F 
R 


CTTGCCAAACTGAGCTTCTG 
GGTGATGGTGCTCCTCTTTC 


62 


150-180 


KF267877 


BA14 


(TC)v 


F 
R 


CAGCAAAGTCCACAGCACA 
TGCGTGCACGTGTGAGT 


62 


140-200 


KF267878 


BA15 


(CA)9 


F 
R 


TCCCATTTTCCTGGTCGT 
TGGCATCCAAGACTCTGCTA 


55 


270-300 


KF267879 


BA16 


(CA)ii 


F 
R 


TCACAGATCCTACAACAGCACC 
AGAAGGAGAACGCAAATACCC 


55 


190 


KF267880 


BA22 


(AG)6 


F 
R 


CCACAGGTTTGTGTTTGTTCTC 
GTCCCATTCCTTTCATATTCCA 


55 


330-360 


KF267881 


BA24 


(GA)3o 


F 
R 


ACAGCGAAGCTCACACACAT 
TCTGTGGAAGAAGAGTGGGAAT 


55 


170-230 


KF267882 


BA25 


(AC)4o 


F 
R 


TGCCCAGTAAATAGGAGAGATTG 
TTATGCTGCTGGAATGGTATTG 


55 


150-180 


KF267883 


BA26 


(CT)25 


F 
R 


AGGT GC AT T GGAAAGGT AT T GA 
GTTTGAGGCTTTGGACATACATC 


65 


360-400 


KF267884 


BA28 


(TG)7 


F 
R 


AGGACTTCCATTTTGTTGGGTA 
GGGTTTTAATTTCGTGACTTGC 


55 


340-400 


KF267885 


BA30 


(CT)33 


F 
R 


CTTGAATGCTGCTGGTAAAGC 
GCATCCTTTTGGACTCAATTTC 


65 


290-370 


KF267886 



Note: = annealing temperature. 



1 min, followed by 30 cycles of denaturation at 94°C for 1 min, 1 min at a specific 
annealing temperature (T^), and a final extension of 72°C for 5 min. The obtained 
products were verified by electrophoresis on 3% agarose gels containing 0.1 mg 
ethidium bromide per milliliter in Ix TBE buffer and genotyped using 6% dena- 
turing polyacrylamide gels dyed with silver nitrate (Creste et al., 2001). We esti- 
mated the allele sizes by comparison to a 10-bp DNA ladder (Invitrogen). 

The amplicons were electrophoretically separated using an ABI 377 auto- 
mated sequencer (AppHed Biosystems) with GeneScan 500 TAMRA marker as 
the size standard (Applied Biosystems). The fragment size and allele identifica- 
tion were determined using GeneMarker version 2.2 software (SoftGenetics, 
State College, Pennsylvania, USA). Cross-species amplifications were evalu- 
ated using five other species from the Rubiaceae family with varying phyloge- 
netic proximity to B. australis: B. mendoncaei K. Schum., B. stipulata (Veil.) 
C. Presl, and Rustia formosa (Cham. & Schltdl.) Klotzsch (Condamineeae 
tribe, Ixoroideae subfamily); Rudgea jasminoides (Cham.) Miill. Arg. (Psycho- 
trieae tribe); and Coussarea accedens Miill. Arg. (Coussareeae tribe, Ru- 
bioideae subfamily) (Bremer and Eriksson, 2009). 



We characterized the preliminary genetic diversity of B. australis populations 
from lowland (20 individuals; 23.3762°S, 45.0806°W, Ubatuba, Sao Paulo, Brazil) 
and upland (20 individuals; 23.3259°S, 45.0710°W, Sao Luis do Paraitinga, Sao 
Paulo, Brazil) Serra do Mar. Descriptive statistics and Hardy-Weinberg equilib- 
rium tests were performed using GenAlEx version 6.5 (Peakall and Smouse, 2006). 
The same nine loci were polymorphic in both populations (Table 2) and for these 
the average number of alleles was 8.6, ranging from five to 13 alleles per locus. The 
observed and expected heterozygosities ranged from 0.333 to 0.900 (average: 
0.629) and 0.564 to 0.900 (average: 0.742), respectively. The fixation index ranged 
from -0.089 to 0.61 1, with an average of 0.147. Four loci in the lower population 
and three loci in the upper population showed significant deviations from Hardy- 
Weinberg equilibrium (P < 0.05). These results indicated some slight excess of 
homozygotes, which might have resulted from mating between relatives and/or the 
inbreeding generated by selfing; B. australis is a self-compatible species (Andrich, 
2008). All 10 loci amplified successfully in the other Bathysa C. Presl species, but 
only four primers performed well for Rustia formosa (Table 3); all primers failed 
for Rudgea jasminoides and Coussarea accedens. 



Table 2. Results of initial primer screening of lower (23.3762°S, 45.0806°W) and upper (23.3259°S, 45.0710°W) Bathysa australis populations. Only 
polymorphic loci are shown. 



Population 


Locus 


A 




Ho 


He 


F 


HWE^ 


Lower 


BA02 


5 


2.417 


0.450 


0.601 


0.232 


ns 


(A^=20) 


BA14 


12 


5.270 


0.842 


0.832 


-0.039 


ns 




BA15 


8 


5.369 


0.850 


0.835 


-0.045 


ns 




BA22 


5 


2.222 


0.500 


0.564 


0.091 


ns 




BA24 


12 


2.462 


0.500 


0.609 


0.158 






BA25 


7 


5.026 


0.353 


0.825 


0.559 






BA26 


11 


8.163 


0.600 


0.900 


0.316 


ns 




BA28 


8 


5.755 


0.900 


0.847 


-0.089 






BA30 


12 


6.968 


0.333 


0.881 


0.611 




Upper 


BA02 


5 


2.606 


0.600 


0.632 


0.026 


ns 


(iV=20) 


BA14 


11 


6.422 


0.647 


0.870 


0.234 






BA15 


8 


3.404 


0.750 


0.724 


-0.062 


ns 




BA22 


6 


2.548 


0.550 


0.623 


0.095 


ns 




BA24 


11 


4.938 


0.800 


0.818 


-0.003 


ns 




BA25 


7 


4.040 


0.750 


0.772 


0.003 


ns 




BA26 


13 


7.407 


0.650 


0.887 


0.249 






BA28 


5 


3.162 


0.450 


0.701 


0.342 






BA30 


10 


4.396 


0.800 


0.792 


-0.036 


ns 



Note: A = number of alleles; A^ = effective number of alleles; F = fixation index; = expected heterozygosity; = observed heterozygosity; HWE = 
Hardy-Weinberg equilibrium tests; A^^ = number of individuals sampled. 

^Significant departures from HWE are indicated at the following levels: *P < 0.05, **P < 0.01, and < 0.001; ns = nonsignificant. 
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Table 3 . Results from the cross-amplification tests using primers designed 
for Bathysa australis. 



Bathysa 



Bathysa 



Rustia 



Coussarea 



Rudgea 



Locus 


stipulata 


mendoncaei 


formosa 


accedens jasminoides 


BA02 


+ 


+ 


+ 




BA14 


+ 


+ 






BA15 


+ 


+ 






BA16 


+ 


+ 






BA22 


+ 


+ 


+ 




BA24 


+ 


+ 






BA25 


+ 


+ 






BA26 


+ 


+ 


+ 




BA28 


+ 


+ 






BA30 


+ 


+ 


+ 





Note: + = successful amplification; - = failed amplification. 

CONCLUSIONS 

The microsatellite markers described here are the first developed 
for the genus Bathysa and will be useful for genetic, ecological, 
and conservation management evaluations. Cross-species amplifi- 
cations suggest that some of these loci may be useful in other spe- 
cies from the tribe Condamineeae (subfamily Ixoroideae). 

LITERATURE CITED 

Andrich, M. 2008. Sistema reprodutivo e polinizagao em duas especies 
arboreas simpatricas de Bathysa (Rubiaceae). Dissertation, Instituto 
de Pesquisas, Jardim Botanico do Rio de Janeiro/Escola Nacional de 
Botanica Tropical, Rio de Janeiro, Brazil. 

BiLLOTTE, N., P. J. R. Lagoda, A. M. RisTERUcci, AND F. C. Baurens. 
1999. Microsatellite-enriched libraries: Applied methodology for 



the development of SSR markers in tropical crops. Fruits 54: 
277-288. 

Bremer, B., and T. Eriksson. 2009. Time tree of Rubiaceae: Phytogeny 
and dating the family, subfamilies, and tribes. International Journal 
of Plant Sciences 170: 766-793. 

Creste, S., a. Tulmann Neto, and A. Figueira. 2001. Detection of single 
sequence repeat polymorphisms in denaturing polyacrylamide se- 
quencing gels by silver staining. Plant Molecular Biology Reporter 
19: 299-306. 

Germano-Filho, p. 1999. Estudos taxonomicos do genero Bathysa 
C.Presl (Rubiaceae, Rondeletieae), no Brasil. Rodriguesia 50: 
49-75. 

Martins, W. S., D. C. S. Lucas, K. F. S. Neves, and D. J. Bertioli. 2009. 

WebSat: A Web software for microsatellite marker development. 

Bioinformatics (Oxford, England) 3: 282-283. 
Myers, N., R. A. Mittermeier, C. G. Mittermeier, G. A. B. Fonseca, 

AND J. Ken. 2000. Biodiversity hotspots for conservation priorities. 

Nature Am-. 853-858. 
Peakall, R., and p. E. Smouse. 2006. GenAlEx 6: Genetic analysis 

in Excel. Population genetic software for teaching and research. 

Molecular Ecology Notes 6: 288-295. 
Pedroni, F. 2001. Aspectos da estrutura e dinamica da comunidade ar- 

borea na Mata Atlantica de planicie e encosta em Picinguaba, Ubatuba, 

SP. Dissertation, Universidade Estadual de Campinas, Campinas, Sao 

Paulo, Brazil. 

Ramos, E., R. B. Torres, R. F. A. Veiga, and C. A. Joly. 2011. Estudo 
do componente arboreo de dois trechos da Floresta Ombrofila 
Densa Submontana em Ubatuba (SP). Biota Neotropica 11(2): 
http://www.biotaneotropica.org.br/vl ln2/en/abstract?inventory+ 
bn024 11022011. 

Temnykh, S., G. DeClerck, A. Lukashova, L. Lipovich, S. Cartinhour, 
AND S. McCoucH. 2001. Computational and experimental analysis 
of microsatellites in rice (Oryza sativa L.): Frequency, length varia- 
tion, transposon associations, and genetic marker potential. Genome 
11: 1441-1452. 



Appendix 1 . Voucher information for species used in this study. 





Voucher specimen 


Collection 


Species 


accession no. 


locality 


Bathysa australis (A. St.-Hil.) Hook. f. ex 


HRCB 60163 


Rio Claro, SP 


K. Schum. 






Bathysa stipulata (Veil.) C. Presl 


HRCB 60107 


Rio Claro, SP 


Bathysa mendoncaei K. Schum. 


HRCB 59786 


Rio Claro, SP 


Rustia formosa (Cham. & Schltdl.) 


HRCB 59785 


Rio Claro, SP 


Klotzsch 






Coussarea accedens Miill. Arg. 


HRCB 59788 


Rio Claro, SP 


Rudgea jasminoides (Cham.) Miill. Arg. 


lAC 49279 


Campinas, SP 



Note: HRCB = Herbario Rioclarense, Universidade Estadual Pauhsta, 
Rio Claro, Sao Paulo, Brazil; lAC = Herbarium of the Instituto Agronomico 
de Campinas, Campinas, Sao Paulo, Brazil; SP = Sao Paulo. 
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